Fluctuation relations without microscopic time reversality: Generalized Green-Kubo 
relation and integral fluctuation theorem for uniformly sheared granular systems 
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We derive the generalized Green-Kubo relation and an integral form of the fluctuation theorem 
that apply to uniformly sheared granular systems in which microscopic time-reversal symmetry 
is broken. The former relation provides an exact representation of nonequilibrium steady-state 
properties as the integral of the time-correlation function describing transient dynamics from an 
initial quiescent towards a final sheared steady state. We also investigate implications of the integral 
fluctuation theorem on the approach towards the steady state and on the possible form of the steady- 
state distribution function in terms of the excess thermodynamic function. 
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Developing statistical mechanics for nonequilibrium 
steady states is one of the most challenging problems 
in physics. Indeed, quite a few studies have been de- 
voted to such a development 0, 0, 0, 0, 0, 0, 0, 0, 0, [HI , 
and among the most remarkable outcomes have been the 
generalized Green-Kubo relation [j| and various forms of 
fluctuation theorems [1, 0, 0, 0] . A suggestive represen- 
tation of the steady-state distribution in terms of the ex- 
cess entropy production has been recently proposed [l(| , 
and its connection to the steady-state thermodynamics 
has been discussed [11] . However, most of these stud- 



ies rest on the use of microscopic time reversality or the 
local detailed balance, and they cannot be applied lit- 
erally to macroscopic dissipative systems like granular 
fluids [13] despite manifest similarities [lj|. In this Let- 
ter, we explore to what extent those outcomes hold also 
for this important class of systems in which microscopic 
time reversality is broken. It is demonstrated that the 
generalized Green-Kubo relation, an integral form of the 
fluctuation theorem, and the steady-state distribution in 
terms of the excess thermodynamic function can be de- 
rived without resorting to microscopic time reversality. 
This will be exemplified for uniformly sheared granular 
systems. Such a formulation will also be useful in studies 
of jammed glassy materials [14] . 

We shall consider a system of N dissipative soft-sphere 
particles of mass m in a volume V subjected to stationary 
shearing characterized by the shear-rate tensor KXfj, — 
jSxxS^y with the shear rate 7. It is assumed that the 
applied shear induces a homogeneous streaming-velocity 
profile k ■ r at position r. Equations of motion describing 
such a system are the SLLOD equations Q, 



hK'f., 

m 

^2 Fij - « • Pi 



(la) 
(lb) 



turn of particle i. We assume pairwise-additivc "smooth" 
contact forces acting only on the normal direction. The 
simplest realistic model for such a force Fy that particle 
j exerts on particle i is the frictional contact model [lit . 

Fy =rijQ(d-rij)[f(d-rij)-j(d-rij)(gij-Tij)]. (2) 

Here d denotes the particle diameter; Q(x) is the Heavi- 
side step function; = r— Tj andgy = (fj— Tj)/m; and 



fij/nj with r*j 



The first term in Eq. 



describes a conservative force representing the elastic re- 
pulsion: typical functional forms are f(x) oc x (linear 
model) and f(x) oc x 3 ^ 2 (Hertzian contact model). The 
second term refers to a nonconservative, dissipative force 
due to inelastic collisions. It is proportional to the rela- 
tive velocity of colliding particles, and breaks the time- 
reversal symmetry of Eq. (llbj) . The amount of energy 
dissipation is characterized by the viscous function "f(x) 
which is typically assumed to be constant or modeled as 
j(x) oc x 1 / 2 . The internal energy of the system is 



H 



(3) 



where the potential energy function V(ry) satisfies 
dV(rij)/drij = —Q(d — rij)f(d — ry). The rate of the 
internal energy change is computed by using Eqs. |T]) as 

H a = -ja xy - 211. (4) 

Here we have used the specific form kx^ — jSxxSfiy and 
the symmetry ax^ = o>a of the stress tensor 



E/ 



(5) 



Here fj and pi refer to the position and peculiar momen- 



7Z is Rayleigh's dissipation function 
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The time evolution of phase variables, i.e., functions 
of the phase-space point T = (r N ,p N ), is determined by 
the Liouville equation [§| 



dt 



A(T) 



d_ 
dT 



A(T)=iCA(T). 



(7) 



The explicit expression for the Liouville operator iC can 
easily be found from Eqs. ([1]). A formal solution to 
Eq. reads A(T(t)) = exp(iCt)A(T). (Hereafter, the 
absence of the argument t implies that associated quan- 
tities are evaluated at t = Q, and the dependence on T 
shall often be dropped for brevity.) On the other hand, 
the Liouville equation for the phase-space distribution 
function f(T,t) is given by Q 



df(T,t) 
dt 



-[i£ + A(T)]f(T,t) = -itff(T,t). (8) 



Here A(T) = (d/dT) ■ T is the phase-space compression 
factor. For the SLLOD equations |T]), one obtains 

A ( r ) = ™ E E Q ( d - r *Md - n 3 ) < 0. (9) 

i j^i 

A formal solution to Eq. © is f(T, t) = cxp(-iCU)f(T). 
In the following, we shall often use the relations [3| 

J dT[e iCt A(T)]B(T) = J dT A(T) [ e - icU B(T)}, (10) 



exp(— iC^t) = exp 



dsA(-a) exp(-iCt). (11) 



Let us consider the following realization of the noncqui- 
librium steady state. The system is first equilibrated at 
the inverse temperature (3 — 1/T (setting Boltzmann's 
constant unity) by turning off the shearing and dissipa- 
tive forces. The distribution function for such a fictitious 
state is given by the canonical one 



/c(T;/3) 



-fiH (T) 



dTe -PH (D_ (12) 



This choice of the fictitious initial state will be justified 
below. At time t = 0, the shearing and dissipative forces 
are turned on, and thereafter the system evolves accord- 
ing to the SLLOD equations |1|. The nonequilibrium 
distribution function for t > is therefore given by 



/(r,t) = e Xp H£tt)/ G (r ; /3). 



(13) 



The steady state is assumed to be reached for t — > oo. 

One obtains from Eqs. (fTl]) and (|13[) the so-called 
Kawasaki representation Q 



f(r,t) 



exp 



ds A(— s) 



-PHo(-t) 



f c (T;P)exp / dsfl(-s) 



(14) 



In the second equality we have used Ho(—t) = Ho(0) — 
L dsHo(-s) and Eq. ((4]), and introduced 

n(T) = ~Pia xy {T) - 2f3K(T) - A(T). (15) 

One easily finds that f2(r) coincides with the dissipation 
function introduced in Ref. Q. From the normalization 
J dT f(T,t) — 1, one obtains a Jarzynski-type equality 
( e / ds n(-s)^p = 1 Q j but a more useful form of such an 
equality shall be derived below. From here on, (• ■ • )p = 
J dT f c (T; 0) ■ ■ ■ refers to the averaging over the initial 
canonical distribution function. 

For the nonequilibrium average (A(t))p defined by 

(A(t)) p = J dTf c (T;0)A(t) = J dT A(0) f(T,t), (16) 



in which the second equality follows from Eq. JT0 
finds using Eq. |T 



(A(t))p = (A(0)eti dsQ ^)t 



one 



(17) 



By differentiating and then integrating this equation with 
respect to time, we obtain 

(A(t)) = (A(0))p + [ ds (A(aM0)b- (18) 
Jo 

The derivation of Eqs. (|17[) and (| 1 8j) in terms of the dis- 
sipation function defined in Eq. (fT5|) is called the dissipa- 
tion theorem [T3| • Although it is implied in Ref. [13] that 
the dissipation theorem is a corollary of the transient fluc- 
tuation theorem resting on microscopic time reversality, 
it holds also in the absence of such a symmetry. 

However, a question arises as to the utility of Eqs. (fT7|) 
and (|18p since they explicitly refer to the fictitious initial 
equilibrium distribution. In this regard, let us prove here 
two important properties concerning lim t ^ 00 (A(t)}/s for 
systems that exhibit mixing [181 ]. First, it follows from 
Eq. (fT8|) for t -> oo 



Jt (A(t)) = (A(t)si(P))f> - (Mt)h(to(o)b = o, (i9) 

meaning that hmt_ >0O (A(t)) / g becomes a constant. Here 
we have used (Cl(Q))p = which can easily be confirmed 
from Eqs. (O, ©, and ©■ Thus, the steady-state aver- 
age (A) ss = lim t _ >00 (A(t))^ is well defined, and is given 
by setting t — > oo in Eq. |[T8|) as 



(A) ss = (A(0)) p + / ds (A(s)fl(O)) . (20) 
Jo 



Second, since 



g -pHa{-t) 
dpi Z(j3) 



} = [{Hoh-H (-t)] 



-PH (-t) 



(21) 
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one obtains from Eqs. (| 14[) and (|16|) for t — > oo 



c9 



dr^(o) [(H ) fj -H Q (-t)]f(r,t) 

= (A(t)) (H Q ) -(A(t)H o )p^O, (22) 



i.e., (A) ss = limt_ >00 (yl(t)) i g is in fact independent 
of the inverse temperature (3 of the fictitious initial 
equilibrium state. Actually, one can prove a stronger 
statement 19]: (A) ss determined from Eq. (|20p co- 
incides with the steady-state average calculated based 
on an arbitrary initial distribution function f(T) via 
(A) ss = ]n^^JdTA(0)exp(-itft)f(F). Thus, (A) ss 
is uniquely specified by the "thermodynamic" parameters 
(N, V, 7) characterizing the steady state, irrespective of 
the initial distribution function. It is therefore most con- 
venient to adopt the aforementioned fictitious state for 
which the canonical distribution function is available. 

Equation l|20p is the generalized Green-Kubo formula 
relating the steady-state average to the time-correlation 
function describing transient dynamics from an initial 
equilibrium towards a final steady state. It is a natu- 
ral extension of the one derived in Rcf. ;|| which takes 
into account effects from inelastic collisions. 

We next derive a generalized Jarzynski-type equality 

Since J* dsfl(-s) = f_ t dsQ(s), the left-hand side can 

be expressed as / df / c (f ; f3) e a ^ dsn(t (s)) . By setting 
f = r(i), there holds, since f (s) = T(t + s), 



(e a fo dsQ (- s )) = J dT(t)f c (T(t);f3)e a f° dsn{T ( s)) . 

(24) 

Using H (T(t)) = H (T) + J*dsH Q (T{s)) and dT(t) = 
e / ds A ( r ( s ))rfr [8] which follows from the conservation 
of the number of ensemble members within a comoving 
phase volume, one derives similarly to Eq. (|14[) 

dT(t)f c (T{t);p)=drf c (r;!3)e-fi d ' a M'». (25) 



Substituting this into Eq. ([24]) yields the equality 

For a = 1, the equality (f2"3")l reduces to the one noted 
below Eq. (fT5"|) . By setting a = in Eq. (f2"3"|) , one obtains 

(e-/o*«*»nW) /3 = i. (26) 

This equality is called the integral fluctuation theorem [9[ 
or the nonequilibrium partition identity [171 ]. originally 
derived as a corollary of the transient fluctuation theo- 
rem that rests on microscopic time reversality. Thus, the 
integral form of the fluctuation theorem holds also in the 
absence of such a symmetry. Using the Jensen inequality, 
one obtains from Eq. 



(27) 



which is referred to as the second- law inequality 17 1 . 

The approach from an initial equilibrium towards a fi- 
nal steady state can be characterized by the change in 
from (Q(0)}p = towards (fi) 88 = -(A) ss > 
[see Eq. ©], which follows from Eq. (|15p since (Hq) ss — 
-j{vxy) ss ~ 2(TZ) SS = [see Eq. (g))]. Notice that the 
inequality (|2T|) implies (Q(t))p > for short t, meaning 
that the system must initially move towards, rather than 
away from, the steady state. This is an analogue of the 
statement originally made for the relaxation to equilib- 
rium from nonequilibrium state [2Cj ■ 

The fact that (D,(t))p > for short t can also be derived 
from Eq. (fTS)) specialized to A = CI, 



(fi(t)) p = / dsmnmp, 



(28) 



since the autocorrelation function must be positive at 
least for short times. If, in addition, the inequality 
(0(4)0(0))^ > holds for all t as, e.g., in an exponential 
decay, one finds from Eq. (2SJ) that {il(t))p < (H(f))p 
for t < t' , and the average dissipation function increases 
monotonically. Since (Vt(t))p is related to the free en 
ergy production rate T{t) via (3J-(t) — {Sl(t))p 



ds (Sl(s))p > for all t 



such a monotonic increase in (Q(t))p implies that the 
steady state can be characterized as the state of maxi- 
mum free energy production rate. However, the inequal- 
ity (f2(t)f2(0))/3 > does not hold for all t in general, and 
the approach towards the steady state is not necessarily 
monotonic. Again, this is an an alog ue of the statement 
for the relaxation to equilibrium [20(. 

Finally, we show that, at least formally, f(T,t) can 
be expressed in terms of the excess quantity. Such a 
representation has been derived exploiting the micro- 
scopic time reversality (lo| , which however is broken here. 
Nevertheless, one can speak of paths evolving forward 
and backward in time whose probabilities are related 
via the integral fluctuation theorem (e^ Jo dsQ (~ s ))p = 
( e -i So d8n ( s )) 01 obtained by setting a = 1/2 in Eq. ([23]) . 
Utilizing such a symmetry, one can derive [22j ] 

/(r,t)// c (r ; /3) = (e*J"o^n(-)) r . /( e -iJ?«tonW) r;tl 

(29) 

where {X) r , t = J dfXf c {f;(3)S(T(t) - T)/f(T,t) de- 
notes the conditioned average. We expand this expres- 
sion using cumulant average (• ■ • )p. t to obtain 

log -||^|y = \ {<e«)r ; o + (e c + x )r ;t } + |(fl- + 0+) 

OO 1 

+ E 2Tm {((©- ) fe )r ; o - (-l) fe ((©+ ) k )h} , (30) 

fc=2 

in which Gf = J*ds[fl(±s) - Q±] with fi± = 

limt^oo(lA) Jo ds (n(±s))p. Notice that {Q^p/fl is the 
excess free energy production (see above). 
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Following Ref. [10j], we introduce a small parameter e 
characterizing the "degree of nonequilibrium" . Express- 
ing all the quantities in natural dimensionless units, we 
shall choose e ~ 7 ~ 7 so that f2 = 0(e) [see Eqs. ([5]), 
©, ©, and dSJ]. Applying Eq. {TSJ| to <fi(±t))0, one 
finds fi_ + = 0(e 3 i) [23]. The factor et comes from 
e iCt _ e i£o* _|_ O(et), where i£o denotes the Liouville op- 
erator at e = 0. Similarly, one finds that the last term in 
Eq. ([30]) is 0(e 3 t). Due to the additional t dependence 
in the second term in Eq. (|30j) . we obtain 

= \ {( 0C - >r;o + (e c + x )r ;f } + 0(6 3 t 2 ), (31) 

which resembles the representation derived in Ref. [loj ]. 
The t dependence is retained in the correction term since, 
e.g., in the non-Newtonian regime where the relaxation 
time towards the steady state can become ~ I/7 ~ 1/e, 
the correction increases to 0(e) on that time scale. Thus, 
the expression (|3"Tj) can be useful only if the steady state 
is reached on the time scale satisfying et <C 1. 

However, such a naive small-e expansion like Eq. (|3ip 
might be questionable. As noted below Eq. (J22J), one 
can choose an arbitrary initial e = state for discussing 
the steady-state properties for e > 0. Then, a question 
arises concerning what e = state the small-e expansion 
should refer to. We would like to leave the discussion on 
the e = versus e — ► state to our future work. 

In this Letter, we derived the generalized Green-Kubo 
relation and the integral fluctuation theorem for systems 
without microscopic time reversality and discussed their 
implications. Our formulation for granular systems will 
also be useful in studies of jammed glassy materials 
In particular, the generalized Green-Kubo relation serves 
as a convenient exact starting point for approximate the- 
oretical treatments, and can easily be combined with 
liquid-state theories as detailed in Ref. [23] ]. 
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